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A method is introduced to optimize excited state trial wavefunctions. Tiie metiiod 
is applied to ground and vibrationally excited states of bosonic van der Waals 
clusters of upto seven particles. Employing optimized trial wave functions with 
three-body correlations, we use correlation function Monte Carlo to estimate the 
corresponding excited state energies. 

1 Introduction 

Solving the Schrodinger equation is a fundamental problem, but also it is a 
timely one in view of the rapid experimental progress made in recent years.EI 
For example, spectra of van der Waals clusters with embedded chromophores 
have been measured and have been used to construct more accurate inter- 
atomic pair potentials and to test proposed three-body potentials. The ability 
to compute ground and excited state properties of clusters is also imnprtant 
for the interpretation of diffraction experiments on small ^He clusters. □ Such 
experiments can be done with transmission diffraction gratings, which have 
become available in recent years. The early experiments showed the presence 
of ^He dimers, and other small clusters, but recently, more detailed physical 
properties such as energies and bond lengths have been extracted from these 
experiments.a In spite of this experimental progress, thus far, no computa- 
tional method exists to answer even the minimalist question of how many 
excited states small ''He clusters have. At a cluster conference in Germany in 
1997, this question made the list of pressing issues. 

Another important issue that can be addressed by studying clusters is the 
magnitude of three-body contributions to the interatomic potential energy. 
Evidence obtained by means of quantum Monte Carlo methods in van der 
Waals complexes such as. Ar„HF and Ar„HF dates back to at least the mid 
nineties.Qtl Hutson et alu studied the role of three-body forces in more detail 
by a recent Discrete Variable Representation (DYp,) study of Ar„HF with 
= 2, 3, 4, and work in this field is still continuing.!!! 

Modifying diffusion Monte Carlo to calculate vibrational states is generally 
considered to be a difficult problem. In special cases, the excited state is the 
lowest energy state of a particular symmetry, and then the standard fixed node 
approximation is applicable. Indeed, this approach has been used to compute 
tunneling splittings for comparison with experiments on some water clusters.I3 
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For electronic excited states, in many cases, approximate wave functions have 
nodal surfaces that are sufficiently accurate that the fixed node diffusion Monte 
Carlo method yields good energies but for many vibrational problems adequate 
nodal surfaces are not available. 

A method with the potentialof addressing the excited state problem is cor- 
relation function Monte Carlo.Ercl A promising new alternative diffusion Monte 
Carlo method for calculating excited states employs the so-called projection 
operator imaginary time spectral evolution (POITSE) approach. Ej An impor- 
tant feature of both of these methods is that approximate knowledge of the 
excited states can be built in from the start, and the statistical accuracy of the 
estimates can be improved dramatically in this way. In fact, only rarely can 
one obtain results of sufficient accuracy without these initial approximations. 
However, generating them can be quite difficult, and this is the problem we 
address in this paper by means of Monte Carlo wave function optimization. 

More specifically, we address the problem of computing energies of vi- 
brationally excited states by means of quantum Monte Carlo methods. We 
propose a method consisting of a combination of linear and non-linear op- 
timization schedules to generate optimized trial functions which are used in 
a correlation function Monte Carlo computation. The method is applied to 
bosonic van der Waals clusters. As do other quantum Monte Carlo methods, 
our approach has the ability to deal with systems with strong anharmonici- 
ties and strong quantum mechanical fluctuations, which cause the failure of 
conventional variational, normal mode, and basis set methods. In contrast to 
the fixed node diffusion Monte Carlo method, the one discussed here does not 
require a priori kaowledge of nodal surfaces. 

The methodtil discussed in this paper relies on the use of optimized trial 
functions for the excited states, and the optimization method is explained in de- 
tail. In applications, once the optimized trial functions have been constructed, 
we use correlation function Monte Carlo to reduce systematically the varia- 
tional bias of the energy estipiates. In principle, the imaginary time spectral 
evolution (POITSE) method,E3 (also see the paper by Whaley in this volume) 
could be used with the optimized excited state wave functions we discuss here. 
It would be interesting to compare the relative merits of these two projection 
methods. 

2 One state 

We consider clusters of atoms of mass fj,, interacting pairwise via a Lennard- 
Jones potential with core radius a and well depth e. In dimensionless form, the 
reduced pair potential can be written as v{r) — r"^^ — 2r^^ and the reduced 
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Hamiltonian a,s H — /2m + V , where P"^ /2m and V are the total reduced 
kinetic and potential energy operators. The only parameter is the dimension- 
less inverse mass-tp"^ = / [icP't, which is proportional to the square of the de 
Boer parameter ,t3 a dimensionless measure of the relative strength of quantum 
fluctuations. 

We use the position representation, and denote by i? the SA^c Cartesian 
coordinates of the TVc atoms in the cluster. Suppose we have a real-valued trial 
function ip{R). Typically, this trial function may have 50-100 parameters and 
it may depend non-linearly on these parameters. First we recall how this wave 
function can be optimized by minimization of the variance of the local energy 
£{R), which is defined by 

HijiR) = £{R)ip{R). (1) 
Following Umrigar et aL,Q one can minimize the variance 

x'^an-in/f), (2) 



which in the position representation can be written as the variance of the local 
energy. Note that is nothing but the square of the uncertainty in the energy, 
so that x^ vanishes for any eigenstate of the Hamiltonian H. 

The minimization of x^ can be done by means of a Monte Carlo procedure 
with the following steps: 

1. Select a sample of configurations Ri, . . . , Rs from the probability density 
V'l (to be defined). 

2. Evaluate: 



B = and B' = I , (3) 



where 





ijg{R) ^g(i?) 



3. Find £ from least-squares solution of 

^'{R^)^£^j{R„), (5) 
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for fj = 1, . . . , s: 

s 

Y,i^{R,)i,'{R,) 
£^^s ■ (6) 

CT = 1 

4. Vary the parameters in the trial function to minimize x^, the normalized 
sum of squared residues defined by the previous step: 

s 2 

E [Vi'(i?.)-Fv3(i?,) 

X' = — -s . (7) 

<T=1 

The least-squares Jiiinimization of is done by means of the Levenberg- 
Marquardt algorithm.tj It should be noted that this wave function optimiza- 
tion algorithm can, in principle be applied to any eigenstate, but the basin of 
attraction of the ground state typically is vastly bigger than that of any other 
eigenstate, so that, the algorithm will practically never produce anything but 
an approximation to the ground state, unless one starts the optimization with 
a carefully designed initial wave function, such as, e.g., described in detail in 
the next section,. 

For the purpose of optimizing only the ground state, the best choice for 
the guiding function tpg, which is used to generate the sample of configurations, 
is the optimized ground state wave function itself. Since this function is only 
known at the end of the optimization, one uses a reasonable initial approxima- 
tion, if available. Otherwise, a few bootstrap iterations may be required. 

For optimization of excited states, one can use a power of the optimized 
ground state trial wave function. We have used a power which is roughly in 
the range from one half to one third. This has the effect of increasing the range 
of configurations sampled with appreciable probability. The goal is to produce 
a sample that has considerable overlap with the all excited states of interest. 

3 Several states 

Next we consider the problem of finding the 'best' linear combination of a 
number of given elementary basis functions /3i, . . . ,/?„. Before we continue, 
we should explain our terminology, since it reflects the procedure that will be 
used. We form linear combinations of the elementary basis functions. These 
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linear combinations depend on any non-linear parameters that appear in the 
elementary basis functions; the linear combinations will be optimized with re- 
spect to the non-linear parameters by means of the general non-linear trial 
function optimization procedure described previously in Section ^ (See the 
discussion after Eq. ([l^) for a more explicit description of the combined opti- 
mization procedure.) Finally, these optimized basis functions which S£Q/e as 
the basis functions in a correlation function Monte Carlo calculation,Q'Ll and 
we shall return to this in more detail later. 

If the 'best' linear combinations of elementary basis functions are defined in 
the sense that for such linear combinations the expectation value of the energy 
is stationary with respect to variation of the linear coefficients, the solution to 
this problem is well known. Being a linear pmblem, the solution requires for 
its implementation traditional linear algebra.Ej The featured matrices consist 
of the matrix of overlap integrals of the elementary basis functions, and the 
matrix of the Hamiltonian sandwiched between them. The trouble, of course, 
is that the required matrix elements can be estimated only by means of Monte 
Carlo methods and the elementary basis functions we employ for the cluster 
problem. 

Stationarity of the energy is equivalent to the least-squares principle that 
is used in the following algorithm. The latter can be used with a very small 
sample of configurations, but in the limit of an infinite sample it produces 
precisely the solution for which the energy is stationary. 

To find the optimal linear coefficients perform the following steps: 



1. Select a sample of configurations . 
ipg (as discussed previously). 



. , Rs from the probability density 



2. Evaluate: 



B 



and 



where 



B' 



\ Pi{Rs) 

( P'liRi) 
\ P'liRs) 

P^{R) 




(8) 



V'g(i?) 



and (3[{R) = 




MR) 



(9) 



(10) 
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3. Find 

from least-squares fit to 



E = 



(11) 



$r^R,) = Y,P,{R,)£,., (12) 

for a = 1, . . . , s and i = 1, . . . ,n. 
4. Find the eigensystem of E and write 

n 

f.,=5:df)^.df , (13) 

k=l 

where the d/'^-' and dj'"'' are components of left and right eigenvectors of 
E with eigenvalue E^. 

This algorithm yields an approximate expression for the eigenstate for energy 
i,^^\R) « V;W(i?) = ^ A(i?)df ^ (14) 

i 

In addition to an approximate eigenstate, this yields.an eigenvalue estimate 
which satisfies the following approximate inequalitjO for the excited state 
energy 

Ek < Ek. (15) 

The inequality holds rigorously in the absence of statistical noise, i.e., if an 
infinite Monte Carlo sample is used, or, if by other means, the quantum me- 
chanical overlap integrals and matrix elements, corresponding to the matrices 
N and H defined in Eq. (p^), are evaluated exactly. 

Before we discuss some technical details of the linear optimization, let us 
summarize the full optimization scheme for excited state k. We iteratively 
minimize the variance of the local energy of the approximate eigenstate given 
in Eq. (p^, where the variance is defined in Eq. (|^). This minimization is 
with respect to the non-linear parameters that appear in the elementary basis 
functions as to be defined exlicitly in Section ^ During the least-squares min- 
imization, for any given choice of non-linear parameters, the linear parameters 

ik') 

dj are defined by the algorithm given in this section, which indeed produces 
the state given by Eq. ( p^ ) required in the definition of the variance in Eq. (^. 

In the ideal case that the basis functions are linear combinations of no 
more than n true eigenfunctions of the Hamiltonian, the previous algorithm 
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yields the true eigenvalues, even for a finite Monte Carlo sample, unless it fails 
altogether for lack of sufficient independent data. 

To find the least-squares solution for E from Eq. (jlj) we write the latter 
in the form 

B' = BE. (16) 
Multiply through from the left by B-'", the transpose of B, and invert to obtain 

E = (B^B)"\b'^B') = N^^H. (17) 

It is simple to verify that indeed this yields the least-squares solution of 
Eqs. (|l2|). Note that the matrix H becomes symmetric only in the limit of an 
infinite sample; symmetrizing it for a finite sample destroys the zero- variance 



propoerty of Eq. 17 



It is is well-known that the solution for E as written in Eq. (|1^) is nu- 
merically unstable .EZl This is a consequence of the fact that the matrix N is ill 
conditioned if the Pk are nearly linearly dependent, which indeed is the case 
for our elementary basis functions. The solution to this problem is to use a 
singular value decomposition to obtain a numerically defined and regularized 
inverse B~^. In terms of the latter, one finds from Eq. (hi 



E = B-iB'. (18) 

More explicitly, one uses a singular value decomposition to write0 

B = USV^, (19) 

where U and V are square orthogonal matrices respectively of order s, the 
sample size, and n, the number of elementary basis functions, while Sr is a 
rectangular s xn matrix with zeroes everywhere except for its leading diagonal 
elements cti > (72 > > 0; r is chosen such that the remaining singular values 
are sufficiently close to zero to be ignored. In our applications, we ignored all 
singular values <7k with < 10"^(Tiedbi, where e^bi is the double precision 
machine accuracy. This seems a reasonable choice, but we have no compelling 
argument to justify it. 

From Eq. ([l^ ) one obtains 

E = V^S^^U^B', (20) 

where U,. is the s x r matrix consisting of the first r columns of U; is the 
n X r matrix likewise obtained from V; and Sr is the r x r upper left corner 
of S. 
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4 Elementary basis functions 

We used elementary basis functions of the general form introduced in Ref. El. 
Rotation and translation symmetry are built into these functions by writing 
them as functions of all interparticle distances. First of all, we introduce 
a scaling function with values that change appreciably only in the range of 
interparticle distances that occur in the cluster configurations with appreciable 
probability. For this purpose we first introduce a piecewise linear function /. 
This function has three parameters: xi < X2 < X3, which define the four linear 
segments of the continuous function /: 



The parameters are determined by the relevant length scales of the system. The 
parameter xi sets the scale for how close two atoms can get with reasonable 
probability; X2 roughly equals the most likely interparticle distance; and X3 is 
the distance at which one expects the onset of the long distance asymptotic 
regime. Possibly, one could drop X2 and use a simpler function consisting of 
three linear segments only. 

The function / has no continuous derivatives and cannot be used directly 
as a scaling function. Instead, we use the generalized Gaussian transform 



with c = 0.1. !_. 

In their most general form the wave functions in Ref.E3 contain five-body 
correlations, but in the work reported here we have only used three-body cor- 
relations and for completeness we shall describe the construction of these func- 
tions explicitly. 

Choose three of the atoms. Suppose they have labels a, (3, and 7 and 
Cartesian coordinates rQ,,r^ and r^. This defines three scaled interatomic 
distances 




1 for X < xi, 

for X — X2, 

1 for X > X3. 



(21) 




(22) 



Define three invariants as sums of powers of these variables 




(23) 




(24) 
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with p = 1,2,3. Clearly, any polynomial in the invariants /i,/2 and I3 is 
symmetric with respect to permutation of the labels a, (3, and 7. A convenient 
property of these variables is that the reverse is also true: any symmetric 
polynomial in the three scaled distances can be written as a polynomial in the 
invariants Ii , I2 and . This makes it simple to parameterize these symmetric 
polynomials. 

In terms of the invariants we define 'minimal polynomials' Si as follow: 
pick a monomial in Ii, I2, and I3 and sum over all possible ways of choosing 
three atoms a, (3, and 7. These polynomials are minimal in the sense that one 
cannot omit any single term without violating the bosonic symmetry. 

In addition to bosonic symmetry, wc impose short and long-distance bound- 
ary conditions. This yields the following form for the elementary basis func- 
tions 



with 




As discussed in detail in Rcf.t2l, the r^T^ term in the exponent and its coefficient 
are chosen so that, when two atoms approach each other, the strongest diver- 
gence in the local energy, i.e. the Lennard- Jones r~^^ divergence, is canceled 
by the divergence in the local kinetic energy. The energy Ek is determined 
self-consistently by iteration; one or two iterations typically suffice. The spe- 
cific form of the decay constant is chosen on the basis of two assumptions. 
The energy is assumed to be proportional to the number of atom pairs in the 
cluster .E2l This is reasonable for small clusters, but for larger ones this should 
probably be modified to reflect the expectation that the energy is proportional 
to the average number of nearest neighbor pairs. The second assumption is 
that if one atom is far away from all others, the wave function can be written as 
the product of an iVc — 1 cluster wave function and an exponentially decaying 
part that carries a fraction of the total energy equal to the number of bonds 
connecting that atom to the others. 

The aj in Eq. ( |25| ) are non-linear variational parameters. Their optimal 
values are re-optimized for each excited state. In principle, one could optimize 
all non-linear parameters, including those that appear in the scaling function 
and the factors that impose the boundary conditions. However, it has been 
our experience that this produces strongly correlated variational parameters 
and results in unstable fits. 
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5 Reduction of variational errors 

The linear and non-linear optimization procedures described above are used to 
generate basis functions for a correlation function Monte Carlo calculation^ 
which increases the statistical accuracy of the energy estimates and reduces the 
systematic errors due to imperfections of the variational functions. The number 
of these basis functions is much smaller than the number of elementary basis 
functions that appear in the linear combinations. The advantage of not using 
all elementary basis functions for correlation function Monte Carlo purposes 
can be understood as follows. 

Suppose that the optimization phase yields states with k = 1, . . . , 

n' < n. Correlation function Monte Carlo in a statistical sense yields the basis 
functions 

|V;«(t))^e-^^|^«). (27) 

As t increases, the spectral weight of undesirable excited states, i.e., states k 
with Ek > En' is decreased. That is desirable, but at the same time all basis 
functions approach the ground state and therefore become more nearly linearly 
dependent. More explicitly, one has Monte Carlo estimates of the following the 
generalization of Eq. ( [l7| ) 

E(i)=N(i)-iH(t), (28) 

with 

iV,,(t) = (^«(t)|^«(0) (29) 

and 

H.,,{t) = {^j;(^Ht)\n\^j^^\t)). (30) 

Again, trouble is caused by an ill-conditioned matrix, which in this case is N(t), 
and increasingly so for increasing values of the projection time t. Obviously, the 
better are the trial states IV'*'*-') and the fewer is their number, the less severe is 
this problem. We should also point out in this context that the singular value 
decomposition cannot be used in this case. The reason is that the analogs of 
the matrices B and B' become too big to store for Monte Carlo samples of the 
size required in the correlation function Monte Carlo runs. 

6 Results 

We computed reduced energies for clusters consisting of various noble gas 
atoms He, Ne, Ar, and Kr, corresponding respectively to the dimensionless 
inverse masses = 9.61 x 10"^^ 7 092 x lO'^, 6.9 635 x lO"'' and 1.9 128 x 
10^**. Quantities with the dimension of energy can be reconstructed from the 
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following values for the corresponding well depths e/ZcbK — 10.22, 35.6, 119.4, 
164.0, 222.3. 

In Table we compare msults obtained with our Monte Carlo method 
with results of Lcitner et aZ.,E3 which were obtained by the discrete variable 
representation (DVR) method. With the exception of the fifth state of Ne, the 
Monte Carlo results agree with or improve the DVR results. In some cases, the 
disagreement can be attributed to lack of convergence of the DVR results. Ea 
Similar effects were observed in prior variational calculations, which employed 
a precursor of the optimization scheme discussed hereJHj The discrepancy for 
the fifth state of Ne may be an illustration of a weakness of the correlation 
function Monte Carlo method, as it is commonly implemented, namely the 
difficulty of estimating the statistical and systematic errors. 

As is well known, in diffusion Monte Carlo computations one has to con- 
tend with time-step errors, due to the short-time approximation that has to be 
used for the imaginary time evolution operator exp —fH.. The same applies to 
the correlation function Monte Carlo method. For all results reported here, we 
repeated the calculations for a range of time steps to verify that the time-stea 
errors were smaller than the statistical and other systematic errors. See Ref. Ilil 
for further details. 

There can be problems both with obtaining reliable estimates of the statis- 
tical errors and with establishing convergence as a function of projection time 
t [c/. Eq. (p8|)]. This is a consequence of the fact that the data for different 
values of the projection time are strongly correlated since they are obtained 
from the same Monte Carlo data. Correlated noise may introduce false trends 
or obscure true ones, a problem that in principle can be solved by performing 
independent runs for different projection times, but that would greatly increase 
the computation time. 

Unreliable statistical error estimates may also come about because the 
correlation fpjption Monte Carlo calculation takes the form of a pure-diffusion 
Monte CarlocSE3 calculation. The algorithm used for the latter features weights 
consisting of a number of fluctuating factors proportional to the projection time 
t. Consequently, as the projection time t increases, the variance of the estii 
mators increases and they acquire a significantly non-Gaussian distribution,E3 
which renders error bars computed in the standard way increasingly mislead- 
ing. Conceivahli^, one could reduce the severity of this effect by using branching 
random walkspO as is dune in standard diffusion Monte Carlo, or by means of 
reptation Monte Carlo.Ea 

In Table || we present results for the energies of the first five levels of Ar 
clusters of sizes four through seven. Our method allows one to go beyond 
seven atom clusters, but, as one can see from Table p|, the statistical errors 
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Table 1 : Comparisson with DVR results of reduced energy of vibrational energy levels 
of noble gas trimers; the estimated errors (which reflect a combination of statistical and 
systematic errors, as explained in the text) are a few units in the least significant decimal. 



k 








MC 


DVR 


MC 


DVR 


1 


-1.719 560 


-1.718 


-2.553 289 43 


-2.553 


2 


-1.222 83 


-1.220 


-2.250 185 5 


-2.250 


3 


-1.142 


-1.138 


-2.126 361 


-2.126 


4 


-1.038 


-1.035 


-1.996 43 


-1.996 


5 


-0.890 


-0.898 


-1.946 7 


-1.947 



increase with system size. To obtain more accurate results for larger clusters 
it would probably be helpful to include higher order correlations in the wave 
function, since the degrees of the polynomials were chosen sufficiently high that 
increasing them further no longer improves the quality of the trial functions. 

Figure |l| contains three energy levels as a function of mass for four particle 
clusters. The harmonic approximation implies that for large masses the energy 
is a linear function of m"^ . We expect the energy to vanish quadratically in 
the vicinity of the dissociation limit. The results are therefore plotted using 
variables that-yield linear dependence both for large masses and for energies 
close to zero.c3 As the energy levels approaches zero, both the optimization 
and the projection methods begin to fail, and correspondingly data points are 
missing. 

In the elementary basis functions, we typically used polynomials of degree 
ten in the prefactors and of degree three in the exponent. The diffusion Monte 
Carlo runs used on the order of a million steps with a time step of a couple 
of tenths. The longer runs typically took a few hours on a four processor SGI 
Origin 200. 

7 Discussion 

We presented a scheme to optimize trial functions to be used in quantum Monte 
Carlo calculations for excited states. We applied this scheme to Lennard- Jones 
clusters, but it can also be applied to systems interacting with more realistic 
pair- and even three-body potentials. Only the elementary basis functions 
would have to be modified, since it should be noted that the contribution in 
the exponent in Eq. (|2^) is specifically tailored to suppress the r~^^ divergence 
of the local energy for the Lennard- Jones potential. 

The method can also be applied straightforwardly to clusters of expcri- 
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Table 2: Reduced energy of vibrational levels Efe of Ar clusters; the estimated errors (which 
reflect a combination of statistical and systematic errors, as explained in the text) are a few 
units in the least significant decimal. 



k 


Ar4 




Are 


Ar7 


1 


-5.118 11 


-7.785 1 


-10.887 9 


-14.191 


2 


-4.785 


-7.567 


-10.561 


-13.969 


3 


-4.674 


-7.501 


-10.51 


-13.80 


4 


-4.530 


-7.39 


-10.46 


-13.74 


5 


-4.39 


-7.36 


-10.35 


-13.71 



mental interest containing chromophores. In that case, more drastic changes 
will have to be made to the elementary basis functions. E.g., the invariants 
defined in Eq. (|2^), which guarantee full bosonic symmetry, are of course no 
longer appropriate. In general, the optimization method as such can be used 
for completely different, even fermionic systems, as long as the elementary 
basis functions are adapted for such applications. 

As far as the projection phase of the computations is concerned, the stan- 
dard implementation of correlation function Monte Carlo in the form of pure- 
diffusion Monte Carlo will become inefficient from systems of increasing size 
and might have to be replaced by a branching algorithm based on branching 
random walks or reptation Monte Carlo.cS In this context also the use of our 
optimization scheme in conjunction with the POITSE method is an interesting 
possibility.EJ 

The most serious problem with the results presented in this paper shows 
up at small masses. Clearly, this is caused by a deficiency in the variational 
freedom of the elementary basis functions employed in the currrent applica- 
tions. The incorporation of four-body correlations, which seems to be all that 
the current basis functions are lacking, is likely to yield more accurate results 
for atoms of small mass. 

Acknowledgments 

This research was supported by the (US) National Science Foundation (NSF) 
through Grant DMR-9725080. It is our pleasure to thank David Freeman and 
Cyrus Umrigar for valuable discussions. 

References 



13 



-0.5 - 



-1.5 



-2 - 



k 


= 1 : 


O 


k 


= 2 : 


+ 


k 


= 3 : 


□ 



o 



-2.5 



Ar 



0.05 



Nc 

0.1 



□ 
+ 

o 



□ 

+ 
o 



0.15 



□ 



□ 



o 



□ 



□ 



o 



o 



o 



o 



0.2 



0.25 



He 
0.3 



0.35 



Kr 



m 2 



Figure 1: —yJ—Ej^ for lowest three vibrational levels (k = 1, 2, 3) of four particle clusters vs 
rn~2. The estimated errors for most energies are smaller, than the plot symbols. Results 
for level fc = 3 become unreliable near He and have not been included. The vertical arrows 
indicate Kr, Ar, Ne, and He; the horizontal arrow indicates the classical value -V6. 
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